Topological surface states in three-dimensional magnetic insulators 
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An electron moving in a magnetically ordered background feels an effective magnetic field that can 
be both stronger and more rapidly varying than typical externally applied fields. One consequence is 
that insulating magnetic materials in three dimensions can have topologically nontrivial properties 
of the effective band structure. For the simplest case of two bands, these "Hopf insulators" are 
characterized by a topological invariant as in quantum Hall states and Z2 topological insulators, but 
instead of a Chern number or parity, the underlying invariant is the Hopf invariant that classifies 
maps from the 3-sphere to the 2-sphere. This paper gives an efficient algorithm to compute whether 
a given magnetic band structure has nontrivial Hopf invariant, a double-exchange-like tight-binding 
model that realizes the nontrivial case, and a numerical study of the surface states of this model. 
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Recent theoretical and experimental work has shown 
that there exist nonmagnetic band insulators in which 
spin-orbit coupling plays a role similar to that of the mag- 
netic field in the integer quantum Hall effect (IQHE) . In 
two dimensions [T], these "topological insulators" have 
robust edge states, observed in HgTe/(Hg,Cd)Te het- 
erostructures |2] , and are predicted to show a spin quan- 
tum Hall effect. The existence of a genuinely three- 
dimensional topological insulator phase |31 S [S] with pro- 
tected surface states, recently observed in Bio.gSbo.i [5], 
is rather surprising because the IQHE does not have a 
fully three-dimensional version, but only layered versions 
of the 2D case. Both 2D and 3D topological insulators are 
nonmagnetic, and in fact unbroken time-reversal invari- 
ance is required for the edge state to remain gapless. The 
edge or surface states of topological insulators and IQHE 
states exist because there are topological invariants that 
distinguish these insulating states from ordinary insula- 
tors, and across a boundary between one of these states 
and an ordinary insulator, the energy gap must close. 

The goal of this paper is to show that there are gen- 
uinely three-dimensional (i.e., not layered) topological in- 
sulating phases of electrons moving in a magnetic back- 
ground. The problem of electrons moving in such a 
background has attracted considerable interest because 
of its relevance to materials in which electrons from 
outer orbitals move in a magnetic environment gener- 
ated by the ordered magnetic moments of core electrons. 
The outer-orbital electrons are frequently described by 
tight-binding models in which the hopping of an elec- 
tron from site i to site j depends on its initial and fi- 
nal spin and the nearby core spin configuration. Several 
materials with layered kagome structures derived from a 
parent pyrochlore lattice, such as Nd2Mo207 [7], have 
been argued via this mechanism to show the quantum 
anomalous Hall effect [8} E] , in which electrons hopping 
in a magnetic background form a two-dimensional inte- 



ger quantum Hall state. We give an explicit cubic-lattice 
model with a nontrivial three-dimensional topological in- 
variant (the Hopf invariant in momentum space) and ex- 
tended surface states and discuss in closing which ma- 
terials might realize the "Hopf insulator" phase. This 
insulator is simpler in some ways than the Z2 topologi- 
cal insulators, since the minimal realization requires only 
two bands (counting spin) rather than four 1^. We ar- 
gue that pyrochlore-lattice compounds with noncoUinear 
magnetic order are realistic candidates for Hopf insula- 
tors, and discuss specific materials in closing. 

First, let us give a general introduction to topological 
band insulators in any dimension. Consider a band insu- 
lator in d dimensions with n filled bands and m empty 
bands. In k space, such a band insulator is described by 
m + n dimensional matrix H(k) which has m positive 
and n negative eigenvalues for any k (assuming Ep = 0). 
Without changing the ground state, we may deform all 
the positive eigenvalues to 1 and all the positive eigen- 
values to —1. Thus H{k) has the form 

H{k) = W{k)I,^^^wHk), (1) 

where is the diagonal matrix with m I's and n — I's 
on the diagonal, and W{k) e SU{m + n). We see that 
for any fixed fc, -ff (k) is a point on the manifiold SU(m + 
n)/Gm.n-, where Gm.n is a subgroup of SU{m + n) that is 
formed by the transformations that leave Im,n invariant. 
We find that G'„,„ = SU{m) x SU{n) x U{1) and 

H{k) G SU{m + n)/SU{m) x SU{n) x Uil) = M^.n- 

So a band insulator is described by a mapping from the 
Brillouin zone T'^ (a d-dimensional torus) to Mm,n- The 
different classes of topological band insulators in d dimen- 
sions are just classes of mappings T'' Mm,n- Time- 
reversal invariance in fermionic systems (or other ad- 
ditional symmetries [10]) imposes additional restrictions 
that can separate topological classes. 
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For two-band insulators, we have Mi.i = S'^. Map- 
pings from — !■ can be non-trivial (since the third 
homotopy group 713(6'^) = Z), and a familiar example to 
physicists is the Hopf map: 

where x'^ = i^Ti^i) ^^nd ktl^ + l-^il^ ~ non- 
trivial Hopf map S*^ implies that there can be non- 
trivial mappings S^. The corresponding two-band 
topological insulators in 3D will be called Hopf insula- 
tors; these insulators must break time-reversal in order 
to avoid Kramers degeneracies at time-reversal-invariant 
points k = — k in the Brillouin zone. Before writing down 
a specific model of a Hopf insulator, we review briefly the 
basic ideas of Chern number and the Hopf invariant and 
give an algorithm to compute the Hopf invariant for any 
two-band model. 

There are two equivalent ways to understand the Chern 
number of a gapped 2D band structure with two bands, 
and both are helpful in understanding the Hopf invari- 
ant in 3D. Such a band structure is given by four real, 
periodic functions of k: the Bloch Hamiltonians are 

H{k) = ai{k)ax + a2{k)ay + aT,{k)<7z + a^{k)l. (2) 

The fourth function a/^{k) plays no role in the topologi- 
cal classification and is omitted in the following, since it 
only shifts the energy levels, but can be experimentally 
important in determining if a material is insulating. The 
gapless condition is ai(fc)^ + a2(fc)^ + a3(fc)^ > 0. Let 
the components of k run from — vr to tt. For each value 
of fc, a direction on the unit sphere is fixed by 

(ai(fc),a2(fc),a3(fc)) 
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Now in two dimensions, the Brillouin zone has the topol- 
ogy of the torus T^, and to classify band structures we 
need to classify maps from to . Because any map 
from the circle to the sphere S"^ can be contracted to a 
point, maps from the torus to S"^ are equivalent topo- 
logically to maps from to S"^ [IT]. These maps are 
classified by an integer "homotopy invariant" that gener- 
alizes the notion of winding number; it counts how many 
times the first sphere wraps the second. Going back to 
the torus, the invariant is 
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where the local current (to be used later in 3D) is 
1 
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The same invariant may be more familiar in terms of the 
ground state spinor \x{kxTky)). While this spinor has a 
t/(l) ambiguity, the "Berry flux" 
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is gauge-invariant and equal to jz- For the three- 
dimensional case, we need to use the classification of 
maps from to 5^ found by Pontryagin [12]. 

The classification of two-band band structures in 3D, 
i.e., maps from to S*^, is complicated because 
includes three 2D tori, and each of these may have a 
nontrivial Chern number. The guess that such a band 
structure is classified by three Chern numbers, which if 
nonzero give layered integer quantum Hall states, plus 
the Hopf invariant is actually incorrect, and in fact maps 
from to do not form a group. Here we will focus 
on the case where the Chern numbers vanish and there 
is no quantum Hall effect. Otherwise the Hopf invariant 
is no longer integer- valued but takes values in the finite 
group Z2.GCD(n,,«^,n,) for Chern numbers [HITS]. A 
geometrical picture of the Hopf invariant is obtained by 
noting that each point on S'^ has a preimage that is a 
circle in T"^, and that the linking number of two such 
circles taken from different points of is just the Hopf 
invariant. For evaluating the invariant, it is easier to use 
the integral, following Wilczek and Zee |14j |24j. 
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Here the gauge field A satisfies the magnetostatic equa- 
tion V X j4 = J. This integral is invariant under a small 
variation Sfi in ^ and hence a homotopy invariant. 

The Hopf invariant is similar to the Z2 invariant or 
Chern parity [TSl [TB] in topological insulators in that its 
standard integral expression uses the gauge-dependent 
quantity A even though the final result Q is gauge- 
invariant when the Chern numbers are zero; alternately, 
one can use an explicitly gauge-invariant but nonlocal 
form |17j . The Abelian Chern-Simons form that appears 
in the integrand gives an invariant that is fully gauge- 
invariant (again, when all Chern numbers are zero) . One 
previous appearance of the Hopf invariant in physics is in 
determining the statistics of solitons in the nonlinear a 
model ; as an example, a magnetic moment field n(r) 
coupled to fermions can have the statistics of its solitons 
changed by the Hopf term generated by integrating out 
"fast" fermions coupled to the field [15] . 

The Hopf invariant can be computed rapidly for any 
two-band problem with zero Chern numbers. We com- 
pute the Hopf invariant numerically by discretizing the 
torus in real space to compute j , then Fourier transform- 
ing and solving the magnetostatic equation in momentum 
space in the discretized version of the gauge 9^^^ = 0. 
This results in an 0{N log N) algorithm for N points 
on the Brillouin zone grid. We studied a particular lo- 
cal tight-binding model as an example of nontrivial Hopf 
invariant, and confirmed that it has the same Hopf in- 
variant and vanishing Chern numbers as the Hopf map 
to S^. Fig. 1 shows the resulting values of Q for grids 
of different sizes for the model 
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FIG. 1: Convergence of the numerical Hopf invariant integral 
([7| to or 1 as number of grid points increases. The parame- 
ter A is the reduction of the hopping elements relative to the 
on-site term in the model Q. The horizontal axis is the linear 
size of the grid {Ntot = N^id)- 



Z| = sin + i ycos{kx) + cos{ky) + cos{kz) — -j , 
= z^a'z, H^va. (8) 

This model is constructed to be continuously connected 
to the Hopf map, but with only low Fourier components 
and hence a local tight-binding representation in real 
space. Its Hopf invariant is +1 (Fig. 1) and Chern num- 
bers are zero. It has an E ^ —E symmetry (since 04 = 
in ([2])) and no states in the band gap —0.25 < E < 0.25. 

As a check, the hopping terms in the model ([s]) can be 
reduced by a factor A while keeping the on-site term 
fixed. This leads to a gapless region for approximately 
0.63 < A < 0.82. To confirm that this gap closing is asso- 
ciated with a phase transition between the ordinary and 
topologically nontrivial states, we computed the Hopf in- 
variant for different values of A (Fig. 1). In the computa- 
tion of midgap surface states below, we also checked that 
there are no midgap states for a boundary between two 
slabs with Ai and A2 either both greater than or both less 
than this gapless range, while the surface state appears 
once one is above and one is below. 

In real space, there are three terms in the tight-binding 
Hamiltonian obtained by computing Fourier components 
of (jsj): an on-site Zeeman interaction —13(72/4; a nearest- 
neighbor spin-dependent hopping of strength 3 A/2; and 
lcngth-2 spin-dependent hoppings of strength A/ 2. These 
specific hoppings and the symmetry of ([s]) are not re- 
quired in order to ensure that the topological invariant 
takes the value 1, and any change to the model that does 
not close the gap will not change the topological proper- 
ties. The topological invariant guarantees that a smooth 
(i.e., adiabatic) boundary between a Hopf insulator such 
as (|8| and vacuum or a trivial insulator will have gapless 
excitations. 
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FIG. 2: (a) The lowest electron and hole bands with sharp 
(001) surfaces in the model Q, for a 16-site-thick slab. The 
Brillouin zone is —n < < n, —n < ky < n. The bulk 
band gap is shown by the dotted lines at i5 = ±0.25, and the 
spectrum is cut off at £ = ±1.5. (b) 3D plot {kx,ky,E) of 
midgap dispersion of (001) surface states near F point. 



We find numerically that a gapless surface state still 
exists even for (at least some) sharp boundaries, as in 
integer quantum Hall and topological insulator states, 
where a sharp boundary between materials of different 
topological invariant generates an edge or surface state. 
We computed numerically the midgap surface states in- 
duced by abrupt boundaries in the (001), (010), and (100) 
directions. In all cases we find gapless, dispersive surface 
states from a single midgap band; there is a single Dirac 
point except for the (001) direction, which has a ring of 
gapless points (Fig. 2). Adding perturbations that break 
particle-hole symmetry shifts the band-crossing energy 
by different amounts for different points on this ring. It 
is not clear that the surface state survives for all sharp 
surfaces, but for the surfaces above, a midgap surface 
state appears precisely when the Hopf invariant changes. 

The extension to multiple bands is more subtle mathe- 
matically for the Hopf insulator than for the IQHE or Z2 
insulator. (Of course, extra bands that do not mix with 
the Hopf-nontrivial bands will not modify consequences 
such as the surface state.) In the IQHE, each band car- 
ries its own integer invariant, with a zero sum rule over 
all the bands [llj . In the Z2 insulator, each Kramers pair 
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carries one Z2 invariant, again with a zero sum rule [1]. 
In these cases, the physical state of an electronic system 
is determined by a sum over the (possibly degenerate) oc- 
cupied bands. However, the generalization of the Hopf in- 
variant, the integration of the nonabelian Chern-Simons 
form, is gauge-dependent once there are degenerate occu- 
pied bands: the integral of the nonabelian Chern-Simons 
form changes by an integer under large gauge transfor- 
mations. 

The multiple-band problem is related to how disor- 
der affects Hopf insulators. For both Chern number and 
Chern parity (the Z2 invariant of a noninteracting dis- 
ordered system |16j). each individual state or Kramers 
pair of a large disordered sample can be assigned one 
invariant, and the state is stable until nontrivial states 
"float" across the Fermi level. If the surface state of 
the Hopf insulator is stable to disorder, then the sur- 
face theory must contain some topological invariant as in 
the 3D %2 insulator's surface [I9j[2n], unlike the 2D case 
where time-reversal symmetry is sufficient |21l I22j . Even 
if the surface state is localized, as likely since there is no 
obvious surface topological invariant, it may still be ob- 
servable by photoemission or tunneling experiments, or 
in clean samples. 

The general requirements for a magnetic material to 
be in the 3D Hopf insulator phase because of spin-orbit 
coupling of itinerant electrons to a magnetic background 
are similar to those required for the 2D quantum anoma- 
lous Hall effect, believed to be relevant to the metal- 
hc frustrated magnet Nd2Mo207 [3 EE]- NoncoUinear 
static magnetic order of core electrons occurs in many 
3D frustrated materials, including the parent pyrochlore 
lattice of the kagome lattice that has primarily been stud- 
ied [H H] (the two-dimensional kagome lattice is obtained 
as one layer from a pyrochlore structure). Although the 
pyrochlore lattice, e.g. in "spin ice" compounds such as 
Dy2Ti207, may lack long-range magnetic order, many 
pyrochlore compounds spin-Peierls distort at low tem- 
perature into a lower-symmetry phase with noncoUinear 
order. 

Hence a realistic microscopic model starts from local 
Hund's rule coupling of outer-orbital electrons to the 
magnetic order of a weakly distorted pyrochlore. Re- 
alization of the specific cubic lattice example discussed 
above seems unlikely in a real solid, although it may be 
possible to simulate this Hamiltonian using the tunable 
spin-dependent hoppings of ultracold atoms in optical 
lattices. An important future direction is to calculate 
the Hopf invariant for effective band structures in frus- 
trated magnetic compounds such as R2M02O7, R a rare 
earth ion. Experimental searches in such compounds for 
surface states using photoemission or tunneling methods 
might reveal the Hopf insulator state directly. CoupHng 
of conduction electrons to complex magnetic order has 
also been discussed in perovskite manganites P5]. 
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